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MESOSCOPIC LINEAR STATISTICS OF WIGNER MATRICES 


A. LODHIA AND N. J. SIMM 


Abstract. We study linear spectral statistics of x Wigner random matrices T-L on mesoscopic 
scales. Under mild assumptions on the matrix entries of H, we prove that after centering and 
normalizing, the trace of the resolvent Tr('H — z)~^ converges to a stationary Gaussian process 
as A —^ oo on scales <C Im( 2 ;) <C 1 and explicitly compute the covariance structure. The 

limit process is related to certain regularizations of fractional Brownian motion and logarithmically 
correlated helds appearing in [34]. Finally, we extend our results to general mesoscopic linear 
statistics and prove that the limiting covariance is given by the A^^^-norm of the test functions. 


1. Introduction 


The goal of this paper is to study the limiting fluctuations as A —>■ oo of the linear spectral 
statistic 


N 

( 1 - 1 ) X^--°{f):=J2f(.dN{E-Xj)) 

i=i 

where Ai,..., Atv are the eigenvalues of an A x A Wigner random matrix T-L. The mesoscopic or 
intermediate scale is defined by the assumption that dpf ^ oo asA^oo, but dAr/A —0 as A — 
oo. Therefore, if / is decaying suitably at oo, only a fraction N/d^ of the total number of eigenvalues 
will contribute in the sum (|l.ll) . In recent years, there has been growing interest in understanding 
the limiting distribution of (HI) on such mesoscopic scales. This interest has stemmed from, e.g., 
the appearance of novel stochastic processes in probability theory [34], conductance fluctuations in 
disordered systems |25l |26| and linear statistics of the zeros of Riemann’s zeta function m , among 
others dsmuEi]. 

Previously, the majority of studies concentrated exclusively on the macroscopic scale where 
djq = 1 and A = 0 in (II.ip . denoted In this case it was proved for many different types 

of random matrix ensembles that, provided / has at least one derivative, the centered random 
variable 


(1.2) := Ar“°(/) - EAr“°(/) 

converges in distribution to the normal law A(0, cr^) as A —)• oo. Furthermore, an explicit formula 
for the limiting variance cr^ was obtained, see [38] l4] . In analogy with classical probability, we refer 
to such results as central limit theorems (CLTs). 

Going to finer scales, the mesoscopie fluctuations of are known to be highly sensitive when 
compared to the macroscopic scale; in fact the CLT must break down if d^ grows too quickly 
|44j . In particular if dTv = A, only a finite number of terms contribute in the sum (HI) and we 
cannot expect a Gaussian limit. The latter case dj^ = N is known as microscopic and will not be 
considered in this article, though see [23] for an extensive review. Before we state our main results, 
let us describe the class of random matrices under consideration. 


Definition 1.1. A Wigner matrix is an A x A Hermitian random matrix W whose entries 
Wij = Wji are centered, independent identically distributed complex random variables satisfy¬ 
ing E|lTjjp = 1 and ERAj- = 0 for all i and j. We assume that the common distribution /x of Wij 
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satisfies the sub-Gaussian decay dfj,(z) < oo for some c > 0. This implies that the higher 

moments are finite, in fact we have ElWjjl'? < {CqY for some C > 0. We denote by 
the normalized Wigner matrix. Finally, in case all the entries Tiij are Gaussian distributed, the 
ensemble is known as the Gaussian Unitary Ensemble (GUE). 

This ensemble was introduced by Wigner who proved that in the limit N ^ oo, the mean 
eigenvalue distribution of the normalized Wigner matrix T-L converges to the semi-circle law. A 
modern version of this result {e.g. Theorem 2.9 in 0) states that this convergence holds weakly 
almost surely, i.e. 

1 1 r2 _ 

(1.3) jv ^ N^qo, a.s. 

'' j=i ^ -^-2 

for all bounded and continuous functions /. In order to state our main Theorem, we need a 
condition on the regularity and decay of the test functions / entering in (|l.ip . For a,/3 > 0, let 
denote the space of all functions with a-Hdlder continuous first derivative such that f{x) 
and f'{x) decay faster than 0{\x\~^~^) as |x| —oo. Finally, recall the notation := 

Theorem 1.2. Let % he a normalized Wigner matrix as in Definition \l.l\ Suppose that d^ = 
where 7 satisfies the condition 0 < 7 < 1/3 and consider test functions /i,...,/m G 
for some a > 0 and fi > 0. Then for a fixed E G (—2,2) in (jl.ip we have the convergence in 
distribution 

( 1 . 4 ) (xr^ifi ),.. .,xr°ifM)) ^ (xih ),.. .,xifM)) 

where (A(/i),... ,X{fM)) is an M-dimensional Gaussian vector with covariance matrix 

(1.5) E(A(/p)A^) = ^ I"" dk\k\ fp{k)Jfi^), l<p,q<M 
and f{k) := (27r)“^/^ f{x) dx. 

This result improves and extends earlier work of Boutet de Monvel and Khorunzhy m who 
proved Theorem II.21 when 0 < 7 < 1/8, M = 1 and f{x) = {x — z)~^ (see also Theorem II.51 belowl. 
Erdos and Knowles proved an analogue of Theorem 11.21 for random band matrices [26] , including a 
bound on the variance of (HH) in the Wigner case with the same condition 0 < 7 < 1/3 |25] . Apart 
from these works, CLTs for (II.ip were also obtained in several other ensembles |13ll51l[Mll21l[T^lllj . 
Although these works extend to scales 0 < 7 < 1, the proofs rely on exact formulas for the 
distribution of the eigenvalues, which are unavailable in the Wigner setting. 

Let us now make some general remarks about Theorem 11.21 On macroscopic scales 7 = E = 0 
the results are different for Wigner matrices since the limiting covariance depends on the fourth 
moment of the matrix entries laiiei. On mesoscopic scales with hxed E G (—2,2), we show that 
this difference vanishes, indicating a particularly strong form of universality for formula (II.5p (see 
also |39|). As with the local regime, the limiting distribution of (jl.ip is universal in the choice 
E G (—2,2) around which ones samples the eigenvalues. Apparently unique to the mesoscopic 
regime, however, is the scale invariance of the limiting Gaussian process: formula (|1.5p is unchanged 
after rescaling the arguments of the test functions by any parameter (see also Section II.ip . Optimal 
conditions on the test functions given in Theorem 11.21 remains a signihcant issue ever since the 
seminal work of Johansson [38|. The latter article suggests that in the macroscopic regime, only 
finiteness of the limiting variance should suffice to conclude asymptotic Gaussianity, see |52| for 
recent progress in this direction. In the mesoscopic regime we believe analogously that optimal 
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conditions for asymptotic Gaussianity of dni) should be that J^\k\\ f {k)\‘^ dk < oo . It is historically 
interesting to remark that (II.5h already appeared in a famous 1963 paper of Dyson and Mehta [23]. 

From a probabilistic viewpoint, the semi-circle law (|1.3I) may be interpreted as the law of large 
numbers for the eigenvalues of a Wigner matrix; this may be considered the first natural step for 
the probabilist. The second natural step is to prove the CLT for the macroscopic fluctuations 
around (|1.3I) . Going now to mesoscopic scales, only the first step has been investigated in detail for 
Wigner matrices, with the corresponding results known as local semi-circle laws, so called because 
they track the convergence to the semi-circle closer to the scale of individual eigenvalues. The local 
semi-circle laws turned out to be a very important tool in proving the long-standing universality 
conjectures for eigenvalue statistics at the microscopic scale im 1211153]. Consequently, a number 
of rehnements of Wigner’s semi-circle law of increasing optimality were obtained in recent years 
[28l[3ll|30l[27]. Such results will play a crucial role in our proof of the CLT at mesosopic scales. 

In order to state the local semi-circle law, it is convenient to work with the resolvent G{z) = 
{T-L — z)~^, lm(z) > 0. Then according to (|1.3I1 . the Stieltjes transform of H 

(1.6) Siv(^) := N-^TvG{z) 


should be close to the Stieltjes transform of the semi-circle: 
(1.7) s(z) := ^y^(x-z)"^^/4’ 


— dx 


The local semi-circle law shows that this convergence remains valid at mesoscopic scales lm(z) = 
0{df^) for 1 <C <C A^. The following is the latest version of this result (in our notation). 


Theorem 1.3 (Cacciapuoti, Maltsev, Schlein, 2014). [171 Theorem 1 (i)] Fixfj > 0 and let z = t+ir] 
with t G M and rj > 0 fixed. Then there are constants Mq, Nq., C, c, cq > 0 such that 

(1.8) P (^|sAr(F; + z/dN) - s{E + z/dN)\ > < {CqY<j"K-i 

for all ^ < fj,\E + t/djsfl < 2 r]/dj\f, > A^o such that ^ > Mq, and q < cq 


To prove Theorem 11.21 we will start by proving it for the special case of the resolvent f{x) = 
{x — z)~^. Indeed, one can interpret sn{E + z/d^) as a random process on the upper-half plane 
El and ask whether, after appropriate centering and normalization, a universal limiting process 
exists. We will show that the function {N/dN){sNiE + z/d^) — Esn{E + z/dN)) converges to the 
r'+ -processes. These are certain analytic-pathed Gaussian processes defined on H. 


1.1. Mesoscopic statistics and regularized fractional Brownian motion with El = 0. 

Fractional Brownian motion is a continuous time Gaussian process Bnit) indexed by a number 
H G (0,1) and having covariance 

(1.9) E(B//(t)B/f(s)) = -I- -\t - sf^) 

where ch is a normalization constant. A generalization of the usual Brownian motion {H = 1/2), 
these processes are characterized by their fundamental properties of stationary increments, scale 
invariance (i.e. Bniat) = a^Bnit)) and Gaussianity. The parameter H is known as the Hurst 
index and describes the raggedness of the resulting stochastic motion, with the limit of vanishingly 
small H to be considered the most irregular (see e.g. Proposition 2.5 in [20]). Although the 
fBm processes were invented by Kolmogorov, they were very widely popularized due to a famous 
work of Mandelbrot and van Ness m and since have appeared prominently across mathematics, 
engineering and finance, among other fields, see [40] for a survey of fractional Gaussian fields. Until 
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recently however, no relation between the fBm processes and random matrix theory was known. 
The latter relation (discovered in [M]) goes via the limit H ^ 0 and the following regularization 


(1.10) <(() := ^ ([«-'• - 11B,(*) + le“‘ - 11B,(*)) 

where Bc(s) := Bi{s) + iB 2 {s) and Bi,B 2 are independent copies of standard Brownian motion. 
One can verify that in the limit rj ^ 0, one recovers precisely the fractional Brownian motion, i.e. 
Bfl\t) = Buit). On the other hand, taking instead the limit H ^ 0 in ll.lOl one obtains a process 
B^\t), about which the following was proved: 

Theorem 1.4 f [34] Fyodorov, Khoruzhenko and Simm). For a GUE random matrix T-Lgve, con¬ 
sider the sequence of stochastic processes 


( 1 . 11 ) 


<"^(t) := log 


det ^^gue - E - 


- log 


det ^^gue - E - ^ 


and W^\t) := VF^^(t) — E,{wji^\t)). On any mesoscopic scales of the form cZat —>■ oo with d^ = 
o(iV/log(A^)) and with fixed r G M, r/ > 0 and E G (—2,2), the process converges weakly in 

L^[a, 6] to B^'^ as N ^ oo. 


In particular, this gives a functional (in L^) version of the CLT of Theorem 11.21 for GUE random 
matrices with fk{x) = log |x — — ir]\ — log |x — ir]\. Either by computing the resulting 

norm ()1.5I) or by computing the covariance of B^'^ as defined in (jl.lOl) . one hnds the logarithmic 
correlations 


( 1 . 12 ) 


E{{B^^\t) - B^^\s)f) = llog + l) • 


(r}) 

Thus Bq ' inherits many of the fundamental properties of fBm, including Gaussianity, stationary 

increments, although now one has the ‘regularized self-similarity’ B^^^\at) = B^\t) (the latter 
following from the scale invariance of the inner product (|1.5I) ). More generally, Gaussian fields 
with logarithmic correlations have received a great deal of recent attention across mathematics and 
physics, see [35l[33] and references therein. The most famous example of such a held is undoubtedly 
the 2D Gaussian Eree Eield (GEE) [50], which has important applications in areas such as quantum 
gravity m, Gaussian multiplicative chaos and Stochastic Loewner Evolution [48] . The GEE is also 
believed to play a central role in random matrix theory. Eor example, similarly to (ll.llh . it has 
appeared in relation to the characteristic polynomial, either explicitly or in what appear 

to be its various one-dimensional slices [37l[34l[56]. More recently it has appeared as the height 
function for the minor processes of random matrices [iniE]. 

Going now to the Stieltjes transform sjsiiz) of Theorem fol the trivial relation ^Re{sAr(r -|- 

ir])} = ^W^i^^(r) suggests that the appropriate limiting object should be related to the derivative 

of Bq^\t). Although such a derivative could obviously be represented by differentiating inside 
the Eourier integral in (ll.lOp . it can also be conveniently represented by a random series. More 
generally, for z G HI and Hurst index H < 1, define the following ‘Cayley’ series 


(1.13) 


r'+ 

^ H 


(z) := 


1 fz-l-i 

71 w 


2H—2 oo 


E 

k=0 


I r(2 — 2H k) f z — i 


T{2-2H)k\ \z + i 

























where are a family of real i.i.d. standard Gaussians. A quick computation with the 

series (jl.l3p shows that it has zero mean and covariance structure 

{1.14) E(r3-(^,)fSG)) = 

with E(r^( 2 :i)r^(z 2 )) = 0. It follows that the processes are stationary on horizontal line 
segments of the complex plane (c./. the stationary increments ( 11.1211 for the integrated version). 
The r''’'-processes were originally introduced by Unterberger [55] in the context of geometric rough 
path theory and stochastic partial differential equations, but since then the relation to random 
matrix theory has apparently gone unnoticed. We will show that Tq"^ (z) is directly related to a 
fundamental object of random matrix theory: the normalized trace of the resolvent. 


Theorem 1.5. Consider the resolvent G{z) = (Ti — z) ^. Under the same assumptions as Theorem 
If.HI the centered and normalized trace 


(1.15) 


Vn{z) 


^ (TV G{E + z/dN) - E TV G{E + z/djv)), Im(z) > 0 

«7V 


converges in the sense of finite-dimensional distributions to rg’^( 2 :) as N ^ oo. That is, for any 
finite set of points zi,... ,zm in the upper half-plane H, we have 


(1.16) {Vn{zi),...,Vn{zm))^{T'^{zi),...,T'+{zm)), N ^ oo. 


Furthermore, the process Vn is tight in the space U(D) of continuous functions defined on a bounded 
N-independent rectangle D C H and Vn converges weakly to Tg"^ in U{D). 


Proof. For the finite-dimensional convergence in (jl.l 6 p , see Section [2j The tightness condition in 
ld{D) follows from Corollary 13.81 □ 


Intuitively, the underlying reason for the covariance structure (11.141) (with E[ = 0) appearing in 
random matrix theory can be traced back to the fundamental relation with the sine-kernel 

- / I sm I TTi r. — To I 1 \ ^ 

(1.17) 


lim lim E(l/Ar(ti -|- ir]i)VN{t 2 -|- ir] 2 )) 
>-0 TV—>-00 


dN=N 


1 sin( 7 r(ti - t 2 )) 

TT 


(G — ^2) 


where (heuristically) going to slightly larger scales djsf = N'^ with 0 < 7 < 1 has the effect of a 
large time separation \ti — 12 \ smoothening out the oscillations in the numerator, thus reproducing 
(I1.14P with H = 0 (see e.g. [I6ll33] for additional heuristics). 

Theorem 01 can now be easily extended to Wigner matrices, starting with the identity 


(1.18) 



Re(l/Ar(i + iv)) dt. 


Next, by the rigidity of Theorem ll.3l we have E|l/ 7 v(i + ^??)P bounded uniformly on compact subsets 
of t and 7 E [<5, 00 ) for fixed d > 0 (see Proposition IB.4p . Then a standard tightness argument (see 
e.g. [36|) combined with (|1.16p allows us to conclude the convergence in distribution as N ^ 00 , 

(1.19) f Ke{V]\f {t + irj)) dt—> f Ke{T'f^{t + ir])) dt = . 

Jo Jo 


This implies that -4- though now in the Wigner case, subject to a more restricted growth 
of the parameter djsf than in Theorem 11.41 In all cases considered here, optimal conditions on the 
growth of disf should be anything asymptotically slower than the microscopic scale, i.e. we expect 
our main results to hold provided only that d]\f = o{N). 
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1.2. Strategy of the proof. Our proof of Theorem ll.5l will follow closely the approach popularised 
by Bai and Silverstein The technique begins by exploiting the independence of the matrix 

entries of T-L to write Tr G{z) as a sum of martingale differences. Then a classical version of 
the martingale CLT implies that only 2 estimates are required in order to conclude asymptotic 
Gaussianity. For the macroscopic regime, this technique was applied successfully to conclude CLTs 
for many random matrix ensembles, though not without significant computations [3 Ha la [Ml urns]. 
The mesoscopic regime is characterized by the situation that Im(z) = 0{d'^) as N ^ oo, which 
is further problematic in that the majority of bounds for resolvents involve powers of Im(z)“^. 
To overcome this we use the rigid control provided by Theorem 11.31 many times, but for technical 
reasons we were not able to avoid obtaining estimates of order Im( 2 ;)“^. Such estimates are 
the source of the restriction on in Theorem 11.21 

To pass Theorem 11.51 onto the general linear statistic dni) of Theorem O we use an exact 
formula (see Lemma lC.il) : 


( 1 . 20 ) 


1 poo poo 

= -Re / / VN{T + i7^)d'i>f{T,rj)dTd7J 

^ Jo J—oo 


where d := ^ and Tj is a certain 2-dimensional extension of /, known as an almost-analytic 

extension US]. Since d'^f is deterministic, we can use our CLT for -\- ir]) to conclude a CLT 
for X™®®°(/). The main problem there is to interchange the distributional convergence for Vy with 
the integrals appearing in (|1.2UI) . To perform such an interchange it will suffice to prove a certain 
tightness condition which will boil down to having sharp control on ]E|V;v(7' + ™ the various 

regimes of r and rj. In the bulk of the Wigner semi-circle with tj/d^ N~^, the optimal bound 
of Theorem 11.31 plays a key role, since earlier estimates involving log(iV) and factors would 
lead to a divergent estimate in the mesoscopic regime. In the regions outside the bulk, or with 
with very small imaginary part rj/d^ <C N~^, we employ the recent variance estimates of |52] (see 
Proposition l3.ll) which have the advantage of holding uniformly in ?/ > 0, but the disadvantage of an 
additional factor appearing in the bound. In this way we are able to remove the assumption of 
very rapid decay, which appears in most studies on the mesoscopic regime mi (min]. In contrast, 
there is no decay requirement in the macroscopic regime and the main important characteristic is 
the regularity of / [52]) while here the decay adds an additional complexity to the problem. It 
remains an interesting open problem to push our CLT closer to optimal conditions on the decay 
and regularity of /, and on the spectral scale d]y. 

The structure of this paper is as follows. In Section [2] we prove the finite-dimensional conver¬ 
gence in Theorem 11.51 on scales 1 <C dAr In Section |3| we extend the obtained results to 

compactly supported functions / € C'c’“(M) and show how to replace the assumption of compact 
support with a suitable decay condition on /. Finally, a brief Appendix is included. 

Acknowledgements: Both authors wish to express thanks to Alice Guionnet for suggesting the 
main techniques used in the paper. The first author in particular wishes to express his gratitude to 
Alice Guionnet, who provided helpful advice and support through the NSF grant 6927980 ^''Random 
Matrices, Free Probability and the enumeration of maps”. The second author wishes to express his 
gratitude to Yan Fyodorov, Anna Maltsev and Jeremie Unterberger for stimulating discussions. N. 
J. Simm was supported on EPSRC grant EP/J002763/1 ^‘Insights into Disordered Landscapes via 
Random Matrix Theory and Statistical Mechanics”. 
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2. Convergence in law of the Stieltjes transform 


The goal of this section is to prove the following: 

Proposition 2.1. Let zi,..., zm be M fixed numbers in the upper half of the complex plane H. 
Under the same assumptions as Theorem 11.51 the function Vn converges in the sense of finite¬ 
dimensional distributions to Tq"^, i.e. we have the convergence in law 

(2.1) {VN{zi),...,VN{zM))^{r'o^{zi),...,r'+{zM)), iv^oo, 
where rQ"^(z) is a Gaussian process on H with covariance C{z\,Z 2 ) defined by 

( 2 . 2 ) ^))2 

and 

(2.3) E(r'+(zi)r'+(z2)) = o 

To prove Proposition 12.11 it is enough to fix a linear combination 

M M ^ 

(2.4) Zm := CpVNizp) = ^ Cp^(Tr(G(£; + Zp/dN)) - ETi{G{E + Zp/dM))) 

p=l p=l 

and to prove that Zm converges in distribution to a Gaussian random variable with the appropriate 
variance. Our starting point is that Zm can be expressed as a sum of martingale differences, to 
which a classical version of the martingale CLT can be applied, see Theorem lA.ll To satisfy the 
conditions of the martingale CLT we shall follow the technique outlined in Chapter 9 of the book 
[6] of Bai and Silverstein, which is equivalent to the work [4]. Our approach is also valid at the 
macroscopic scales considered in [4] and we feel gives a somewhat more accessible proof in this case. 


2.1. Method of martingales and some preliminaries. Here we outline the martingale method 
and provide the notation used in the remainder of this Section. Let denote the conditional 
expectation with respect to the u-algebra generated by the upper-left k x k corner of the Wigner 
matrix W. Then we have the martingale decomposition 

N 

(2.5) Zm = Xk,N 

k=l 


where 

( 2 . 6 ) 


M 


Xk,N ■= (Efc - Efc_i) ^ Cp-j—Ti:{G{E + Zp/dN)) 


p=i 


Therefore, to prove Proposition 12.11 it will suffice to check the following two conditions: 

(1) The Lindeberg condition: for all e > 0, we have 

N 

(2-7) ^E(|Xfc^Ar|^l|Xt, jv|>e) 0’ N ^ QO. 

k=l 

(2) Conditional variance: we have the convergence in probability 

N M 

(2.8) Efc_i[|Xfc^jv|^] —)• ClCmC {zi , Zm), N ^ oo, 

k=l l,m=l 

N M 

(2.9) ClCmC{zi,Zm), N ^ QO, 


k=l 


l,m=l 
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where C{zi, Zm) = {i{Ti — Tm + iir]i + f?m))) ^ denotes the covariance in Proposition 12.11 


Before we proceed with the proof of these conditions, we provide some of the relevant notation. 


Important notation: Until now the complex numbers Zp were independent of N. For notational 
convenience and in the remainder of this Section only, we will now allow the implicit A^-dependence 

( 2 . 10 ) 

dN 

As before, the sequence d]\f ^ oo as N ^ oo with d^/N 0 and Tp,rip are fixed real numbers 
with rjp 7 ^ 0. We fix FI E (—2 + 5,2 — 5) strictly inside the support of the limiting semi-circle for 
some small <5 > 0. 

Let Hk be the — 1 x — 1 Wigner matrix obtained by erasing the row and column from 
H. We denote by Gk{z) = {Hk — z)~^ the corresponding resolvent. The following formula, a 
consequence of the Schur complement formula from Linear Algebra, will play an important role; 


( 2 . 11 ) 


Tr(G(z))-Tr(Gfc(z)) 


l + hlGk{z)^hk 
T~lkk ■ 2 ' hp,Gi;(^z)hk 


where hk is the column of % with the entry removed. 

Recall the following standard notation for convergence of random variables in L^. For a sequence 
of random variables we write Xjv = Olp{u{N)) to mean there exists a constant c such 

that E|X 7 v|^ < cu{N) for all N large enough. We will repeatedly use the standard fact that if Xj^ 
converges to X in probability and Yat converges to zero in p > 1, then Xi\f + converges to 
X in probability. 

We start with the proof of the Lindeberg condition (j2.7h which follows from the following stronger 
result (due to the trivial inequality | jvl>e — ^‘^~^\^k,N\^ with 5 > 2): 


Lemma 2.2 (Lyapunov). For all mesoscopic scales 1 <C djq <C with e > 0, there is an integer 
5 > 2 such that 


( 2 . 12 ) 


N M 

^ E (Efc - Efc_i) ^ Cp^ Tr G{z) 

p=i 


k=l 


S 

—)■ 0 , 


N ^ QO. 


Proof. By the triangle inequality it suffices to verify the claim when M = 1, ci = 1. By definition 
of Efc we have (E^ - Ek-i)dfi^TiG{z) = (E^ - Ek-i)Zk,N where Zk,N ■= dff^{TicG{z) - TrGkiz)). 
Then Schur’s complement formula implies 


(2.13) 

1 1 + hlGk{zi)‘^hk 

^kk - zi - hlGk(zi)hk 

(2.14) 

= :w(l + (^i) + Tr(Gk(zi}^}Gkk(zi} 

dN 

where we made 

use of the identity for the diagonal elements of the resolvent 

(2.15) 

Crkk(,^l) t / \ 

T~^kk hi^Gk\Z\)hk 

and defined 


(2.16) 

4"(zi) := h\Gk{zi)^hk - N-^ Tr(Gfc(zi)") 








By the conditional Jensen inequality, we have Hence it is sufficient to prove 

N 



k=l 


The limit (j2.17p follows from standard concentration inequalities applied to the variables {z)^ 


Gkk{z) and dp^N ^ Tr G{z)‘^. In particular, Lemmas IB . 1 [ IB . 2l and 12 . 71 show that for any fixed q > 0, 


we have the estimates 


(2.18) 

(2.19) 

( 2 . 20 ) 


d],^6^/{z)=OL.{{dN/Ny/^), 

Gkk{z) = 

d]^^N-^TTG{zf < Ol. (max{d-^ (d^/iV)n) 


Then applying Cauchy-Schwarz and choosing J > 0 large enough, we obtain (12.1711 . 


□ 


Remark 2.3. In the macroscopic regime dp^ = 1, one can argue similarly that Zk^N = (1 + 
s'{z)){—z — s(z))“^ + with high probability. The leading term in this asymptotic is 

deterministic and does not contribute to (E^ — Efc_i)Zfc^ 7 V) while the error term is small enough to 
imply (I2.17p . 

We now proceed to the remaining and most challenging part of the proof of Proposition 12.11 
which is to verify condition ()2.8p . Before we proceed, it’s worth noting that both and \Xk\‘^ are 
finite linear combinations of terms of the form 


(Efc - Efc_i)^TrG( 2 ;i) x (E^ - Efc_i)^TrG(z 2 ) 


( 2 . 21 ) 


and so it suffices to prove the convergence for a single mixed term in the linear combination. Setting 


Yk{z) := (Efc -Efc_i)d^^TrG(z), 


( 2 . 22 ) 


our essential goal in the remainder of this section will be to prove that we have the convergence in 
probability 


N 


1 


(2.23) Gn{zi,Z2) := '^V.k-i\Yk{zi)Yk{z2)] 


N ^ oo. 


(f(Tl -T2 + i{'ni +??2)))^’ 


k=l 


In what follows, the proof of (12.2311 is divided into 3 main subsections: in section 12.21 we rewrite 
Gm{zi.,Z 2 ) in a form suitable for the computation of asymptotics, then in section [ 23 ] the main 


asymptotic results are obtained and finally in section [2~il they are used to prove ()2.23ll . 

2.2. Simplifying the covariance kernel. Our first Proposition shows that Gpf{z\, Z 2 ) can be 
approximated in the following way 

Proposition 2.4. In terms of the variables (12.1611 . define the covariance kernel 




Then we have 
(2.25) 
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Proof. As in the proof of the Lindeberg condition, we start with Schur’s complement formula which 
implies that 

1 l + hlGk{z)^hk 


(2.26) 


Yk{z) = (Efc 


- z - h\Gk{z)hk 

Rewriting Yk{z) via the small terms (I2.16P and expanding, we obtain the exact identity 


Yk{z) = (Efc -Efc_i) 


(2.27) 
where 

(2.28) ek,N{z) := (E^ - Efc_i) 


d 1 -Hkk-S^/iz) 


+ ^k,N{z) 


dzdNZ + N ^Tr(Gfc(z)) 

1 {nkk-8f{z)fGkk{z) 1 6f{nkk-sf{z)) 


dN {z + iV-1 Ti{Gk{zW dN {z + AT-I Ti{Gk{z))) i ' 


This identity is implicit in the work [3] (see Section 4.1 in [3]), but we provide the derivation in the 
Appendix, Lemma IB. 3 1 Then as in the proof of (j2.17l) . we see that ek^Niz) = O^i (dN/N) uniformly 
in k. Similarly, by Lemma lB.il we can replace {z + N~^ Tr(Gfc(z)))“^ with — 5 ( 2 ;), costing an error 
of the same order. Therefore, we have Yk{z) = Yk{z) + 0]^i{dM/N) where 


(2.29) 


Yk{z) = -(Efc - Efc_i)^^s(z)(iRfcfe - 4'W) 


Using properties of the conditional expectation, we compute that 
N 1 92 

(2.30) '^¥.k-i[Yk{zi)Yk{z2)] = s{zi)s{z2) + Cn{zi,Z2) 

k=l Af 1 2 

Then the covariance Gn{zi,Z 2 ) can be estimated as 


(2.31) 


N 


N 


Cn{zi,Z2) = Cn{zi,Z2) + ^Efc_i(Yfc(2;i) - Yk(zi))Ykiz2) + '^'Ek-i{Yk{z2) - Yk{z2))Yk{zi) 


k=l 


k=l 


N 


+ ^Ek-i{Yk{z,) - Yk{zi)){Yk{z2) - Yk{z2)) + 0^1 ( 1 / 4 ) 


k=l 


where we used that s'{zi)s'{z 2 ) are uniformly bounded for any fixed E € (—2, 2). By our estimates 
for Yk — Yk, the last term above is Oii{d'j^/N). For the middle terms, we apply Cauchy-Schwarz 
(twice) to obtain 

N N 


E 


k=l 


^Ek-iYk{z,){Yk{z2) - Yk{z2)) < E^ |n(zi)||yfc(z2) - ^(^ 2 )! 


k=l 


(2.32) 


< 


N 


N 


Y,nyk{zi) -Yk{zi)\^Y.^\^k{z2)\^ 

\ k=l 


k=l 


, ^E\Ykizi) - Ykizi)\^^ECNiz2,z^) + d)^4'4)|2 

\ k=l 

The first term in the product above is 0{^Jd^qjN). In the remainder of this section, it will become 
clear that Cn{zi.,Z 2 ) is bounded in L^, see Proposition 12.171 □ 
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Remark 2.5. Before we proceed further, note that the approximate covariance kernel Cm{zi,Z 2 ) 
in (12.241) is naturally expressed in terms of the auxiliary kernel 

N 

(2.33) Kn{zi,Z2) := 

k=l 


Then by Cauchy’s integral formula and analyticity, we can write the derivatives in (I2.24p as 

(2.34) ^n(zi,Z 2) :=(f dui i du)2-, -- —KNiuJi,uj2) 

d% { 2 my Is Is (zi - ujiy{z2 - uj2y 


where is a small circle with center z and radius l/{2dN )\/(ti — r 2 )^ + (r/i — 7 / 2 )^- This ensures 
that for fixed Ti,T 2 ,rii,r ]2 and N large enough, 5^^ and Sz 2 are disjoint sets. In the degenerate case 
that zi = Z 2 , it is enough to use the Cauchy integral formula with a single circle S^i ■ 

If we obtain uniform estimates on Kjsf of the form Ki\[ = Ki\f + Oii{u{N)), then the error in 
approximating Cn{zi, Z 2 ) is of the same order in N: 


E 


I 


1 


(2.35) 


< 


d% (2vri)2 
1 I 
d% 47r2 


dioi ® duj2 


s{u}l)s{uJ2) 


Js: 


'^2 


dui ® 

Js 


dL02- 


{zi - ujiy{z2 - uj2y 

|s(a;i)||s(a;2)| 


\Kn-Kn\ 


22 


Zi - wiP|z2 - a;2P 


u{N)\ 


IGvr^ (n - T2y + im - r]2y 

The conclusion of this remark is that it will be sufficient just to understand the convergence in 
of the kernel K]\f{zi, Z 2 )- 


Our first Lemma in this direction rewrites K]sf{zi, Z 2 ) in terms of the matrix elements of the 
resolvent Gk{z) := {Tik — z)~^. We will frequently make use of the shorthand notation := 
Gk{zp), p = 1, 2 to emphasize the dependence on the variables zi and Z 2 - 

Lemma 2.6. The covariance kernel Ki\f{zi, Z 2 ) satisfies the exact identity 

N 

(2.36) KNizi,Z2)=N-^Y.^k-i Y. 

k=l i<kj<k 

N 

(2.37) + N-^ Y ^k-iY^k{Gi^^)A{G^k^)iiPik 

k=l i<k 


where /3ik is expressed in terms of the fourth moments 

(2.38) Afc := E(|ITifc|2 - 1)2 
Proof. By definition we have 

(2.39) Efc(4^(zi)) = iV-i ^ Ek{G^^^)ijW;fWjk-N-^Y^k{G^k^)jj 

i<k,j<k j<k 

. .uTTu/.. _L 


(2.40) 


N-^ Y Ek{GyyijW,kWjk + N-^Y^k-iiGyyjy\w,k\^ -1) 

j<k 


i<kj<k 
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Multiplying out the resulting terms, we obtain 

(2.41) Efc_i[Efc(4'=(z))Efc(4'^(z))] 

(2.42) =N-^Ek-i Y1 MG^k^)ijmG^k\M^w,kw;;;:wgk 

i<kj<k 

p<k,q<k 

(2.43) +N-^Ek-i ^kiG^k^hW;^Wjk'^Ek{G^^^)j,{\Wjk\^ 

i<k,j<k j<k 

(2.44) +N-^Ek-i Y1 ^k{Gf\jWkWjkY.Ek{G^^^)jj{\W,k\^ 

i<kj<k j<k 


(2.45) 


+ N-^Ek-i Y. mG^^hEk{G^‘^^)jj{\Wjk\^-l){mk\^ 

i<k,j<k 


1 ) 

1 ) 

1 ) 


It is clear that the sums (I2.43h and (j2.44jl are identically zero, since the vector {Wik}ij^k consists 
of centered independent random variables satisfying E|Vbjfcp = 1. Similarly, the first summation 
(I2.42P will be zero unless i = q and p = j, this gives the hrst sum on the right-hand side of (12.3611 . 
The last term ()2.45l) will be zero unless i = j, which gives the second sum in (|2.37p . □ 


We now proceed with the estimation of the two sums in (j2.36p and (j2.37p . To do this we need 
precise estimates on the resolvent matrix elements appearing in the sums. 

Lemma 2.7 (Bound on the resolvent). Let z = E + as in (j2.inp and consider an off-diagonal 
resolvent matrix element {Gk{z))pq with p ^ q. Then for any positive integer s we have positive 
constants c, G such that 

(2.46) E|(Gfc(z))p,r <(Cs)“(|^) ' 

for all A:,p, g, E, r E R, r/ > 0 and > 0. For the diagonal matrix elements, the same result holds 
but with (Gk(z))pp — s(z) in place of {Gk{z))pq. 

Proof. This follows from Lemma 5.3 in m, see also previous works on the local semi-circle laws 

[2911281I3Q1EI]. □ 

The second sum (I2.37P over diagonal elements of the resolvent can now be dispensed with im¬ 
mediately: 


Lemma 2.8. Assume that 
(2.47) 


sup E(|ITifc|^ - 1)^ < oo 
N,i,k>0 


Then for all points zi and Z 2 with non-zero imaginary part and on all mesoscopic scales 1 <C ^ 
N, we have the convergence in L^: 

N 


(2.48) 


1 

d% dzidz2 


s{zi)s(z2)N ^ EE Ek{G^^)nEk{Gf^)iAk 


0 , 


N ^ oo 


k=li<k 

Proof. By Lemma 12.71 we have 

N N 

(2.49) N-^ EE Ek{Gi^^)iiEk{G^^^)iiPik = N-hizi)siz2) EE Pik + d]\f/N) 

k=li<k k=l i<k 
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where we assumed that the fourth moment of each matrix entry Wit is finite. Since the left-hand 
side of (I2.49P is analytic, we can use the strategy outlined in Remark 12.51 Hence after inserting 
(I2.49P into the left-hand side of (j2.48l) . the Oii{^dN/N) bound remains of the same order while 
the leading term is of order since s and s' are analytic and uniformly bounded provided 
E G (—2 + 6,2 — 6). This completes the proof of the Lemma. □ 

Remark 2.9. This shows that in the mesoscopic regime, the process V]sf{z) is insensitive to the 
value of the fourth cumulants of the matrix elements, provided they are finite. This is in contrast 
to the regime of global fluctuations where the fourth cumulant is known to appear explicitly in the 
limiting covariance formula, see e.g. This suggests that the Gaussian fluctuations obtained 

in the mesoscopic regime are more universal than in the global regime. 

The remaining challenge now is to compute the first sum in (I2.36p . For that we shall show that 
to leading order on the scales 1 <C <C the sum 

(2.50) Si^k{zi,Z2) ^ ^ ^k{Gk{zi))qp^k{Gkiz2))pq 

p<k,q<k 

satisfies a self-consistent equation. The proof relies on heavy computations. 

2.3. Calculation of the sum Si^kizi, Z 2 )- The goal of this section will be to prove the following 

Proposition 2.10. Consider the quantity Si kizi, Z 2 ) in (I2.50p . We have the estimate 

(2.51) 

N — k N — k 

ziSi^k{zi,Z 2 ) = -s{zi)Si^k{zi,Z 2 ) - j;^s{z 2 ) - j;^s{z 2 )Si^kizi,Z 2 ) + OLi{{d%N~^)^G'j 

where the O^i bound is uniform in k = 1 ,..., N . 

In what follows, the proof of this Proposition will be divided into a number of Lemmas which 
eventually culminate in Lemma (j2.16p . We start with the Green function perturbation identity 

(2.52) ziGkizi) = -In-1 + ^ Wijeie'^Gk{zi) 

i^k 

where the vector e* has entries {ei)j = Sij for j < k and (ei)j = dij-i. It turns out to be helpful 
to separate out the correlations between Wij and Gk by introducing the matrix Gkij defined as the 
resolvent of the matrix Hk with entries {i,j) and {j,i) replaced with 0. Simple algebra shows this 
is a perturbation of the original resolvent: 

(2.53) a[" - G<;> = . 

where Cij = 1 if f / j and cu = 1/2. We now insert ()2.53l) directly into (|2.52l) leading to the 
following expansion. 

Lemma 2.11. The matrix elements of the resolvent (G'^^^)pq satisfy 

^iiG^\q = -6pq + iV-1/2 ^ ITp,(G« ).-? - s(^i)(l - 3/(2iV))(G«)p, 

(2.54) '^pj I ^pj I (iGl,pj)jj — s{zi)){Gl^)pq — N '^Cpj{\Wpj\ — l)s(zi)(G^ ^)pq 

j^k j^k 
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Now inserting (12.5411 into the definition (12.501) . we obtain the decomposition 

(2.55) ziSi{zi, Z2){k.) = + 5i^2 + + 5i^4 + 5i^5 + 5i^6 

where 

(2.56) 5i,i = -iV-i ^Efc(Gf))pp 

p<k 

(2.57) S,,2 = N-^/^ Y. 

p<k,q<k 

(2.58) 5i,3 = -s(^i)(1-3/(2A^))5i 

(2.59) 5i,4 = -N-^ Y - sizi)){G^Yp<iMG^k'^)qp 

p<k,q<k 

(2.60) 5i,5 = -iV-2 ^ Cp,Efc(|hrp,p-l)s(zi)(GW)pgEfc(G(')),p 

i¥^k 

p<k,q<k 

(2.61) Si,6 = -lV-" E 

p<k,q<k 

To estimate these sums we will apply Lemma [2.71 repeatedly, together with the following Lemma 
which allows us to bound the matrix elements of the perturbation {Gkij{z))pq. 

Lemma 2.12 (Bound on the perturbation). Provided that W has finite moments of order 4s, the 
same bound (j2.46p of Lemma [T71 holds with {Gkij{z))pq in place of iGk{z))pq. 

Proof. Using the perturbation identity (|2.53l) we get 

(2.62) {Gkij)pq = {Gk)pq + CijWijN ^^^{Gkij)piiGk)jq + CijWjiN ^^‘^{Gkij)pj{Gk)iq 
Iterating (I2.62p once more and applying Minkowski’s inequality shows that 

\iGkij)pq\^ < Cs\iGk)pq\^+Cs\W,j\^N-^G (\{G k)p^{G k) jq\^ + {{G k)pj (G kU^ 

+ Cs\Wij\'^^N ^ {\{Gkij)pi{Gk)ji{Gk)jq\^ + \{Gkij)pj{Gk)ii{Gk)jq\^ 

p\iGkij)piiGk)jj{Gk)iq\ P \{G kij)pj{Gk)ij {G k)iq\ ) 

for some constant Cs depending only on s. From the deterministic bound | (Gfcjj)pq| < HGfcijH <dN/rj 
and ()2.46p . we find by Cauchy-Schwarz that 

(2.63) E\{Gkij)pqf < CsE\{Gk)pq\^ + + 0((djv/(??iV))*) 

where we used that E|Wjj|^® is bounded. Note that the obtained error term is smaller than 
{dN/{r]N)YG found in (12.4611 . This completes the proof of the Lemma. □ 

As is suggested by the structure of the terms (I2.56ll - (|2.61l) . in what follows we will encounter 
many sums of the following generic form 

(2.64) <^(VF) n 

Ci.C2v,Cr me/ 

for some index set I and a complex valued function (f. By simply counting the number of occurences 
of we obtain 
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Lemma 2.13 (Trivial bound). Suppose that E|i?i)(lT)p is uniformly bounded. Then 

(2.65) E|0| < ciV^(d7v/iV) 1^1/2 

with the same result holding if any of the factors in (I2.64p are replaced with the perturbation 
for any a,b £ I. 

Proof. This follows by repeatedly applying the Cauchy-Schwarz inequality in conjunction with 
Lemmas 12.71 and 12.121 □ 

We start by giving some first estimates on the terms to defined in (I2.56I) - (I2.61I) . 

Lemma 2.14 (Estimates of 5iq, 5i^3, 5i_4 and We have the following bounds, holding 

uniformly in = 1 ,..., A^: 

(2.66) 5i,i = -{N- k)s{z2) + Oil((d^iV-i)i/2) 

(2.67) 5i, 3 = -s(zi)5i + Oil {dN/N) 

(2.68) E|5i^4| < c^Jd^N~^ 

(2.69) E|5p6| < C'\Jd%N-^ 

Proof. These estimates are straightforward consequences of the trivial bound of Lemma 12.131 □ 


With 5 i^ 5 we have to take a bit more care because the trivial bound actually produces a bound 
of order d^ which is divergent. To fix this we have to exploit independence and the fact that 
E|Wp,f = 1. 

Lemma 2.15 (Estimate for 5i,5). We have the following bound, holding uniformly in A: = 1,..., A^: 

(2.70) E|5i,5| < c^Jd%/N 


Proof. First denote 

(2.71) Rp:=^{Gi\,Ek{G^^% 

q<k 


and note that the usual estimates imply that E|72pp = 0{d'j^). Now by Cauchy-Schwarz we can 
estimate 5i^5 as 


(2.72) 


E|5i, 5| < Ar-2^ 

p<k 

= n-^Y. 

p<k 


EliJppE ^ CpjiCpj2{\Wpji\^ - l){\Wpj2\^ 
jlj^k,j2^k 

\ 1/2 

jj^k ) 



1/2 


□ 


To complete the proof of Proposition 12.101 the main task is to control 5i_2. The problem is that 

(2) (2) (2) 
the matrix G\. still depends on Wpj, so in the following we will replace G), with G\^y Remarkably, 

the leading order contribution will come from the error in making this replacement. 

Lemma 2.16. We have the following bound, holding uniformly in k = 1,... ,N: 

(2.73) 5i,2 = -{N - k)N-hiz2)Si,k{zi,Z2) + OLi{{d%N-y/^) 
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(‘ 2 \ ^o\ 

Proof. We start by replacing (G^ ^)qp with {G)^p-)qp in the definition ()2.57l) of 5i^2 and denote the 
modified sum by ^ 1,2 • We will show that 81^2 converges to 0 in . We have 


(2.74) 


e|£i.2 |" = iV-’E 

P2,Ql,92 


X Efc {G^kp 2 j 2 )i2g2®"fc (^1^*272)92^2 


where here and unless otherwise stated, all indices in the summation run from 1 to A; — 1. To make 
and Gkp 2 j 2 independent we use a similar perturbation formula to remove the matrix element 
from Gkp 2 j 2 - Therefore, we define resolvent of the matrix IT with the 

column and row erased and with entries and Wj^p^ replaced with 0. It is easy 

to show that we have a similar identity 


(2-76) - 

Then if we replace the last two factors in (I2.74p with (I2.75P the main term has identically zero 
expectation unless both pi = ji and p 2 = j 2 , which we assume for the moment does not hold. The 
higher order terms in (|2.75n give rise to an error 


(2.77) A = N-y^E ^ Wp,j,Wp 2 MG^l^^)nMG^^^^Xp^ 

jlj2,Pl 

P2,gi,g2 

(2.78) X 

and three further error terms which do not differ in any important way from x4 above. To estimate 
A, we now force independence in lTp 2 j 2 by replacing with Again the main term 

obtained by this replacement has expectation 0 because EITp 2 j 2 = 0. The error terms in (I2.75h give 
rise to sums of the form 


(2.79) 71' = JV-"E ir,.,.»--p,,,Et(G«f«),,,.Ei(»'„,JG®™“),.,^^ 

Jlj2,Pl 

P 2 ,gi,g 2 

( 2 . 80 ) X 

We call such a term maximally expanded because we can no longer exploit independence of the 
different factors to reduce the size of the sum (none of the factors have zero expectation at this 
point). See also [32] for related methods. By the prescription (I2.65p . we have 

(2.81) A' < cN-^N^{dNN-^f < cd%N-^ 


It is clear that all error terms resulting from the replacement (I2.75P give the same bounds after 
employing this procedure. 

Now consider the diagonal terms pi = p 2 and ji = j 2 contributing to ^ 1 , 2 : 

(2.82) JV-=>E |M-',,|7Ei(GW)„.Et(G«),„E,(0®.),„EaG«^)„, 

< cN-^N'^idNN-^f < cd%/N 


(2.83) 
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We conclude that E|£li^ 2 p < cd%/N. Now again using (I2.53h we have 

(2.84) 5i, 2 = iV-3/2 ^ + Op{d%N-^) 

j<k 

p<k,q<k 

(2.85) =-JV-^ c«»rtEl(0W)^,Et((0™)„(0f)*) 

j<k 

p<k,q<k 

(2.86) -N-2 ^ 

j<k 

p<k^q<k 

By (|2.65|) . we see that (|2.85|) is d^j^N~^). The final term (j2.86p (let us denote it 5i,2) can be 

written 

(2.87) S,,^ = -N-h(z,) Cp5Ei(GW)^,Et(Gi5,|,.)„. 

j,P,<l<k 

( 2 . 88 ) - N-h{z2) 

j,p<k q<k 

(2.89) -iV-2 ^ Cp,|Wp,fE,(GW),,Efc(G(2))^^.((G®)pp-^ 

j,p,q<k 

To estimate the term (|2.88p , first replace G). with Gj^^j and note that by (|2.65p the error in making 
this replacement is Op{disfN~^^‘^). We denote the remaining sum .Fi^ 2 . We have 

(2.90) E|J-i,2|2 = N-M^2)f Y. 

j2 ,P2,g2 

(^■^^) ^ ^k{G\.p^j^)j2q2^k{Gl^p^j^)q2j2 

This sum has a similar structure to that seen in the computation of E|Ti^ 2 p. Applying exactly 
the same procedure shows that E|.Ti^ 2 p < cd^/iV^. The term p2.89p is djqN~^) as follows 

directly from the generic bound (|2.65p . In (|2.87l) all terms in the sum where p = j will only 
contribute 0^i(dArN“^) so can be neglected, which justifies us setting Cpj = 1 in (I2.87p . 

Finally, note that (I2.87P has a similar form to Si^k except with the perturbed resolvent elements 

(G^kpj)jq- By the trivial bound p2.65l) they can be exchanged with the original ones {G^j}^)jq at a cost 
0/^1 ((dAriV“^)). Then the summation over p just gives a prefactor (A^ — k), leading to (I2.73h . □ 

2.4. Proof of Proposition 12.11 With the most significant challenges dealt with in the previous 
subsection, our aim now is to complete the proof of Proposition 12.11 by solving relation (I2.5ip and 
computing the limit N ^ oo. In other words, we finally prove the conditional variance formula 
P2.8p which is enough to verify Proposition 12.II 

Proposition 2.17. Consider the approximate covariance kernel Gn{zi,Z 2 ) given by (|2.24l) . Then, 
in terms of the notation (|2.10l) . we have the estimate 

Gn{zi,Z^) = + OLl(l/div) + OlI (log(d^7)y^) 
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(2.92) 






































and 


(2.93) 


CNizi, Z 2 ) = OLi{l/dN) + 0^1 


(log{dN)^Jd%/N^ 


Proof. Recall that C'Ar(zi, ^ 2 ) can be expressed in terms of the auxiliary kernel Kj\f{zi, Z 2 ) in (I2.36D : 


(2.94) 


Cn{zi,Z 2 ) 


32 “^- ^s{zi)s{z2)Kn{zi, Z2) 

dj^ dzidz2 


In turn, K]\f{zi,Z 2 ) was expressed in terms of the sum 51 ^^( 21 , 22 ) defined in (I2.50p and computed 
in Proposition I2.1UI Let us denote the error term in that Proposition by Then we have 


(2.95) 


5l,fc(2l,22) 


N-k s{zi)s{z2) _ EN,k _ 

N 1 - E^s{zi)s{z2) 1 - ^ 5 ( 21 ) 5 ( 22 ) 


where E]^^k = Oii{y d^/N) uniformly in k. Then by definition (I2.50p the kernel K]\f{zi,Z 2 ) can 
be written as 


(2.96) Kn{zi,Z2 ) 


1 

N 


E 


N-k 5 ( 21 ) 5 ( 22 ) 

N 1-N^sizi)s{z2) 



k=l 


_ EN,k _ 

1 - ^S{ZI)S{Z2) 


+ OLl{dN/N) 


These sums are close to Riemann integrals, but before we calculate them we must take into account 
a subtle feature of the mesoscopic regime: when the points 21 and 22 have opposing signs in their 
imaginary parts, there is a singularity in the denominators of (j2.96p . due to the asymptotic formula 

(2.97) 1 5 ( 21 ) 5 ( 22 )-^^ EOid^) 

If the signs of the imaginary part are the same, there is no singularity and one finds that the limiting 
covariance is identically zero, as in (|2.93p . Now let us control the errors in (|2.96p . assuming there 
is a singularity. For 0 < u < 1, let ip{u) := |1 — tt5(2i)5(^)|“^. Then standard results about the 
Riemann integral show that the error term in ()2.96p is bounded in by 


(2.98) fjfE- g < sfEE [I Hu) du + 




N 


where ||' 0 ||tv is the total variational norm of the function if . Since ^5(21)5(22) has non-zero 
imaginary part, exists and is Riemann integrable. Then the total variational norm can be 

written 

(2.99) ||i/,||tv = [ \ij'{u)\du 

Jo 

and a simple calculation taking into account the asymptotics (I2.97P shows that HV'IItv = 0{dN) as 
N 00 , while Jq du = 0{log{d]s[)) as N —)• 00 . 

The error from approximating the first term in (j2.96p can be estimated in the same way, and we 
get 


( 2 . 100 ) 
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Inserting (I2.10nh into the definition of Z 2 ) and bearing in mind Remark 12.51 we have 

1 

(2.101) Cn(zi,^) := ,2 a a 

dzidz2 

(2.102) = (log(<i«)v/rfp^) 

(2.103) = (i(T, _ + i(„ + nt))? ^ (>oe(‘iA')v^4/") 

where in the last line we applied formula (|2.97l) and used that 

(2.104) s'{z^)s\-^) = + 0{djl). 

This completes the proof of Proposition 12.171 Therefore we have finally verified the conditions (12.8p 
and ()2.9p . hence completing the proof of Proposition 12.11 □ 

3. Tightness and linear statistics 

The goal of this section is to extend the CLT obtained in the previous section to a wider range of 
test functions, focusing on the distributional convergence of m- Our approach will be to combine 
the Helffer-Sjostrand formula (11.2011 with the main result of the last Section, telling us that the 
hnite-dimensional distributions of converge to Tg"^. If we can interchange the integrals in (11.201) 
with this distributional convergence then the CLT for mesoscopic linear statistics (II. ip would be 
proved. 

In practise we will need to prove a certain tightness criteria in order to justify this interchange, 
involving uniform estimates on the second moment K\Vn{t + ir])\'^. We will handle those estimates 
mainly with Theorem \L 2 i but this Theorem only applies in the bulk region with r] > dj^ /N. For 
smaller rj we need the following variance bound, due to Sosoe and Wong (in our notation). 

Proposition 3.1. Let e > 0. Then for 0 < r] < 1, and \E + T/d]sf\ < 5 we have that there is a 
universal constant C > 0 such that: 

(3.1) E\Vn{t + zr?)p < Cd%r]~‘^~^, 

Proof. See Proposition 4.1 in [52] . □ 

The reader should compare this result with that obtained from Theorem 11.31 which gives an 
improved bound E| l/)v(t + *i?)P < Cr]~‘^ (Proposition IB.4]l but in a more restricted region. 

To prove the CLT for (II.ip . we consider two situations. Firstly we consider the case that / G 
C'c’“(M), the Holder space of compactly supported functions / such that /' is Holder continuous with 
exponent a > 0. The hypothesis of compact support will allow us to avoid complications coming 
from the edges of the spectrum and will serve as a warm up, illustrating our general approach. 

In the last subsection we consider the more challenging case where the hypothesis of compact 
support is replaced with a more general decay condition f{x) and f'{x). In particular we will prove 
our main result. Theorem [L21 

3.1. Compactly supported functions with Holder continnons first derivative. To obtain 
the CLT we will apply the Helffer-Sjostrand formula (|1.20p with 

(3-2) 4'/(t, v) = (fit) + i{f{t + v)- f{t)))Jiv) 

where J{r]) is a smooth function of compact support, equal to 1 in a neighbourhood of r/ = 0 and 
equal to 0 if r/ > 1, see Lemma lC.il It follows that 

(3.3) d^f{t, rj) = (fit) - f'{t + rj)+ i{f'{t + r/) - /'(t))) J(r?) + (ifit) - (fit + rj) - f{t)))f{rj) 
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and if we assume that f' is Holder continuous with exponent 0 < a < 1 then we also have the 
bound \d'i>f{t,r])\ < cr]^ for some constant c > 0 for 77 small. Also '^f{t,r]) satisfies = f{t), 

and that d'iff{t,rj) is compactly supported whenever / is. 

Theorem 3.2. Suppose that satisfies the condition 1 <C ^Ar <C and consider compactly 

supported test functions /i,..., /m whose first derivatives are Holder continuous for some exponent 
a > 0. Then for any fixed E G (—2, 2) in (11.11) we have the convergence in distribution 

(3.4) (^r^°(/i),..., ^r^°(/M)) ^ (Xi,..., xm) 

where (Xi,... ,Xm) is an M-dimensional Gaussian vector with covariance matrix 

1 r°° . — 

(3.5) E(XpX,) = —y \k\fp{k)ffik)dk, l<p,q<M 


Proof. Since / has compact support, we may write 


(3.6) 


1 


1 /•T+ 


Xr^°(/) = -Re 


TT 



0 Jt- 


+ iT])d^ f{T, rf) dr dij 


for some fixed r_ and r+. First note that the region 0 < r] < d^/N can be neglected, since 


(3.7) 

(3.8) 


I 


dn/N j'T+ 


IE(|VAr(T + Z77)|)|(9^'/(r,?7)| dr, dtj < cd% 


f 


dN/N i'T+ 


^OL—l—e 


dr dr] 


= c(r+ - T-)N''{dN/Ny 


The latter goes to zero after choosing 0 < e < (1 — 7 ) 0 , where dj\f = N'^ with 0 < 7 < 1. Therefore, 
it suffices to study the convergence in distribution of 


(3.9) 


X{VnJ) = -Re 

TT 


1 rT+ 


0 Jt- 


Vn{t + ir])xN{'n)d'i'fiT, rf) dr dr] 


where XNiv) is an indicator function, equal to 1 on the region d^/N < 7 < 1 and 0 otherwise. We 
wish to show that XiV^, f) converges in distribution to X(rQ'^, /). Since we proved in the previous 
section that the finite-dimensional distributions of VAr(r -|- irf) converge pointwise to rQ'^(r -t- irf), 
we can appeal to Theorem lA.21 Let denote the set of functions (^{z, tc) : H x C ^ C of the form 
4>{z,w) = d^g{z)w where 7 is a function in the class stated by Theorem 13.21 and M is as in the 
Theorem IA.2I 

Let D be the domain [0,1] x [r_, r+]. Then Theorem [X2] guarantees the convergence in distribu¬ 
tion {X{VnJi),...,X{VnJm)) ^ (X(ro+, /i),..., X(ro+, /m)) provided we check the following 
tightness conditions: 

(3.10) inf limsupP I — Re [ \V]\r{Tir])xN{r])9'^ dr dr] > e\ =0 

ScH,A(S)<oo N^oo Jd\B J 

(3.11) lim limsupP (—Re / {\Vn{t ir])xNi'n)9'^ ~ dr dr] > e] =0, 

K^oo Jd ) 

First we prove (I3.in[) . By Markov’s inequality it suffices to check that 

(3.12) inf limsup / E{\Vn{t ir])\)xN{v)\d^ dr dr] = 0 

BcD,X{B)<oo N^oo JD\B 

Now since r is fixed, the real part appearing in the denominator of the resolvent is bounded away 
from the edges (due to the compact support of /). Also the imaginary part appearing in the 
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denominator is no less than \/N. Hence we can apply Proposition IB.4I and obtain ¥.\V]sf{T + ir])\ < 
cri~^. By the assumptions on /, r]~^\d'I> f{T,r])\ is integrable on D\B, and so ()3.12p is bounded by 

(3.13) inf / cri~^\d"^ f(T,r])\dT dn = 0 

BcD,X(B)<ooJd\B 

where the last equality follows from the dominated convergence theorem. To check (13.1111 we 
proceed similarly, noting that for any 6 > 0 , we have 

(3.14) {\Vn{t + i'n)xN{v)d'^f{T,v)\ - K)^ < ^ + ^^)Xiv(??)94'/(r, 77 )]^+'^ 

Then (13.lip is bounded by 


lim K ^ lim sup 


(3.15) 


K^oo 


N^oo 


< lim K ^ 

K^OO 


/ cirj 
JD 


E( I Vat (r + 777) I )xAr id) \d^ fir, r])\ drj 

Id'if fir, 77)1^"'"'^ drdr] = 0 


D 

- 1-6 


where in (|3.15ll . we choose 5 so that 0 < <5 < q:/(1 — a) so that the above integral is finite by 
the behavior |947y(r, 77 )! < cr]~°‘ for small 77 . This completes the proof that XiVi^^f) converges in 
distribution to the random variable 


(3.16) X(rg+,/) = -Re / rQ+(r + 777 ) 9 T/(r + 777 ) dr d 77 , 

^ JD 

Since the integral of a Gaussian process is Gaussian, it just remains to compute the covariance and 
verify it is well-defined, see the next Lemma. □ 


Lemma 3.3. Consider the random functional defined by: 

(3.17) X(T'^, f) =—Re [ {t + ir])d'if f{T + irj) dr dr]. 

JB. 

Then the covariance EX[T'^^, fi)X[T'^^, f 2 ) is given by (11.51) . provided that /i and /2 are in (^^’“(M) 
for some a > 0, and |/j|, |//| are 0(|x|“*^^'''^)) for some /3 > 0. 

Proof. Recall that the process T^J" (r-|- 777 ) for 77 > 0 appearing in (|3.17l) has the covariance structure 

(3.18) 


E(r'o^ in + ir7i)r'o^ in + iv 2 )) = 


iiin - t 2 + 7(771 772)))^ 


while E(rQ’^(ri -|- 777 i)rQ’^(r2 -|- 7772)) = 0 . We compute 

(3.19) E(x(r'+,/i),x(r[+,/ 2 )) 


(3.20) 


(3.21) 


=^E 

27r^ 


-|- Re 


[ dzidz 2 { Re [V'^^izi)r'f^iz 2 )d'i/f^izi)d^f 2 iz 2 )] 

JHxH I 


T’+iz,)d^f,iz,)K^in)d^hin) 


we proceed to check that the expectation is finite and that we can exchange expectation with 
integrals by utilizing Fubini’s theorem. It suffices to verify that 


(3.22) 


/ 

JHxH 


dzi dz2E\T'^f izi)T'ff iz2)\\d^fiizi)\\d^f.2in)\ < 00, 


by using Cauchy-Schwarz, the above integral is bounded by 


(3.23) 


/ 

JHxH 


dzi dz2rii S2 fiin)\\d^f2iz2)\ < 00, 
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since by our assumptions on fi, we have 

(3.24) \fl(t + r^)-m\ < min(C'i7?“,C2|tl-(i+^)) < 


for any a G (0,1) and c > 0 a constant independent of t and r/, giving us the bound 
(3.25) \d^f.{t,r])\ < 


taking a < /3/(l + /3) gives a bound that is integrable (recall that ^^f^(t,r]) is supported on 
r] G [0,1]). Hence, the covariance of X{r'^,f) is finite, which implies, in particular that X(rQ'^,/) 
is finite almost-surely — yielding that the process is well-defined as required. 

After taking the expectation inside the integral, the first term (|3.20l) vanishes identically. In the 
second term (I3.2ip . we write the covariance in the Sobolev form 


(3.26) 


IEro+(r-Fz7?)r[+(T-Fz7?) = ^ [ dA:|/c|?v)^i(fc)f;^(/c) 

Jr 


where rT-^^jj^(x) = {x — ti — ir]i)~^. We now interchange the integration over k with the integration 
over H X H. To justify this, note that 7 g 7^(^) = —and due to the conditions on 
/i) / 2 , we have the bound: 

(3.27) 


/HxH 


dzidz2 / (zi)| 15^2(^2)1 <C 


1 rl 


drji dr ]2 


0 JO 


(hi + mY 


< oo. 


After interchanging these integrals, the integration over El x El factorises as a product. Now we 
would like to interchange the Fourier transform with the integral over El, for which it suffices to 
bound the following 


(3.28) 


/*1 poo poo 

/ dr]i / dx dri 

JO j —OO J —oo 


54'/(ri,?7i) 


X — Ti — ir]i 


= / C?hl||5r;i * Vlll 
Jo 


where g-rjiiri) = \Ti+ir]i\ ^ and h^i(ri) = |9'kj(Ti, r/i)|. To bound the norm of the convolution, 
we apply Young’s inequality Wgrj^ * hr^-^ ||i < Wgrj^ ||p||^??i ||g with q = 1 — 6 and p = 1 + 6/{l — 5) with 
(5 > 0. A simple computation shows that \\grj^ ||p < cpi while for sufficiently small <5, ||q is 

bounded uniformly in gi due to the integrability assumptions on / and its derivatives. This shows 
that so that p.28p is finite. 

After performing all such interchanges of integration, we finally obtain 

(3.29) E(X(r'+,/i)X(r'+,/ 2 )) 

(3.30) =— / \k\dk / dxe~^^^— / dzi - 94'/-,(zi) 

^ ^ 87rJ_J' y_oo ttU x-zi 

^ i 7 I I 

(3.31) X / dxe~^^^— / dz 2 - f^(z 2 ) + c.c. 

7-00 TT Jh X- Z2 

Now the inner integrals over H can be evaluated by Lemma 1C.11 There is a caveat however, 
Lemma 1C.11 requires the function / to be compactly supported. We remedy this by taking our 
function / and multiplying it by a cutoff function (/>„ = (p{x/n) where (p{x) is 1 on [—1,1] and 
vanishes outside [—2,2], we let fn = 4>nf- By Lemma [CT] we have the identity 

2 poo POO 2 

(3.32) - dgi j dri - —d^fY'^i,gi) = fn{x)+iH[fn]{x). 

TT Jo 7-00 X-Tl-ipi 

It is well known that 77 is a bounded operator from L^(R) to itself, therefore if we take the limit 
as n ^ oo on both sides of (I3.32p (and note that /n —t / pointwise everywhere and in L^(M)) 
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have 


1 


(3.33) 


lim 

n^oo 




dri - 

-oo x-Ti- ir]i 


—d^fn{n,Vi) = f{x) + iH[f]{x), 


we check that the limit interchanges with integral by noting that 


(3-34) \d^uit,r])\ <2\f^{t) - + r])\\J{r])\ +2(|/(t)| + |/(t + ? 7 )|)| 

(3.35) - f^{t + 77)1 <\(j)nit + V)\\f{t + V) - fit)\ + \(Pn{t + ri)- <(i)ll/(OI 

(3.36) + \(j)n{t + 77)||/'(t + r/) - f{t)\ + \f{t)\\(j)n{t + 7?) - </>n(t)|, 


and that ^n(i) = is inifinitely differentiable and therefore in C^’"(M) for any fe G N, a € 

(0,1). So we may bound for large t, 154' f^{t, r])\ by Cr]^^\ for some constant C and any 

a G ( 0 , 1 ), it follows by the dominated convergence theorem that we may interchange limit with 
integral to obtain 

1 /•OO roo 1 

(3.37) - dr]i dri -^54'/(ri, 771 ) =/(x) + 7 i 7 [/](x). 

TT Jo 7-00 X - Ti - 7771 

The Fourier transform of the Hilbert transform is given by 

(3.38) H[f]{k) = -isgn{k)f{k). 


Hence, inserting (I3.37|) into (I3.3ip and applying (j3.38p yields the limiting covariance structure 

1 


(3.39) 


IE(X(/i)X(/ 2 )) = ^ / fiA:|fe|/i(A:)/ 2 (/c)|l- 7 sgn(/c)pd/c + c.c. 


1 


' —OO 
/•OO 


= —J dk\k\h{k)f2{k) 


□ 


3.2. Functions supported on the real line. The main goal of this subsection is to remove 
the assumption of compact support from the functions / in Theorem 13.21 subject to the following 
decay condition on / and /': for some /3 > 0 and |x| large enough, /(x) and f'{x) are 0{\x\~^~^). 
This will complete the proof of the corresponding statement in our main Theorem 11.21 We begin 
by approximating / by a compactly supported function whose support grows at a rate 0{dM) as 
N ^ 00 . The support of the test function can now extend over the edges of the spectrum, so that 
the resolvent bounds of Proposition IB. 41 cannot be applied. The goal of this subsection is to prove 
the following 

Theorem 3.4. Let d^ = N'^ with 0 < 7 < 1/3. If for some a > 0 and P > 0, we have f G C'^’“(M) 
where f{x) and f'{x) decay faster than for large |x|, then the random variable X™®®°(/) 

converges in distribution to a Gaussian random variable with variance given by (1331) of Theorem 
\3.^ Moreover, the multidimensional version stated in Theorem \3.2\ continues to hold for functions 
of this class. 


We will prove this Theorem by means of the following Lemma and two Propositions. 


Lemma 3.5. Let 0Ar(x) denote a smooth cutoff function equal to 1 in when \E + x/d^l < 2 and 
equal to 0 when \E + x/d]\f \ > 4. Let fN{x) := f{x)(pN{x). Then X^®®°(/) = -T™®®°(/jv) + 0 ^ 2 ( 1 ). 


Proof. We follow the same technique here as in Section 4 of |52j . Note that X“®®°((1 — 4>N)f) is 
only non-zero when Ai <—4orAAr>4, so that 


(3.40) 



((1 - T'A')/) (dN^i) 



<P(Ai >4)+P(A7 v < -4), 
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these quantities are bounded by e by [301 Lemma 7.2], Thus 


(3.41) 

(3.42) 

(3.43) 


r^\\f\\L°° 

E | X ^-°(( 1-<^^)/)|2 = 2 / xP 

Jo 


< 




^P 


y^(i - (i>N)f{dNXi 

i 

y~l(l - 4>N)f{dNh] 

i 

\ —>• 0 ; 


> X dx, 


> 0 


□ 


Now we apply the Helffer-Sjostrand formula ()1.20l) to the function /jv, obtaining 

(3.44) =-Re [ [ Vn{t + ir])xN{T)d^ dr dr] 

Jo Jr 

where xn{t) is the indicator function of the region \E + r/dj^l <4. Repeating the derivation of 
EH we see that 

rdff/N /■ _ 

(3.45) / / E{\VN{T + iri)\xN{r))\d^f^{T,T])\dTdr]^0 

Jo Jr 

on all scales of the form dw = N"^ with 0 < 7 < 1. To apply Theorem IA.21 we replace the N- 
dependent /at with /, noting that 54'/^^ - ST/ = ST/^_/, where / 7 v(r) - /(r) = /(r)(^Ar(r) - 1) 
is supported on the region {r/disf ~ d^\ >4. 


Proposition 3.6. 

(3.46) 


ITe have 



TT 



Vn{t + ir])d'i’f^_f{T, r]) dr dr] = 0 ^ 1 ( 1 ). 


We postpone the proof of this Proposition until the end of this subsection, where it will follow 
from a more general argument. Thus Xj 


^meso 


(3.47) 


Iif,VN) = -Re 

TT 



(/) = L(/, Vat) + 0 ^ 1 ( 1 ) where 
Vn{t + ir])XNi'n)XN{T)d'iff{T, rf) drdr]. 


To show that /(/, Vn) converges to /(/, Tq"^), thus completing the proof of Theorem 13.41 it remains 
to check the following tightness result. 


Proposition 3.7. Consider the domain D = [0,1] x M. We have the following estimates: 

(3.48) inf limsup / E{\VN{t + ir])\xN{r])XN{r))\d^ f{t,r])\dr] dt = 0 , 

BcH,A(B)<oo N^oo JD\B 

(3.49) lim K~^ lim / E(| VAr(t + ir/)|^+^X 7 v(r/)xAr(r))|ST/(t, 7y)|^+'^ dry dt = 0. 

K^oo N^oo J£) 

Proof. First we give some bounds on |ST/(t,ry)|. By our assumptions on /, we have 

(3.50) \f{t + ry) - f{t)\ < min(Ci7y“, < cry^^r 

for any a E (0,1) and c > 0 a constant independent of t and ry. Hence by construction we have the 
bound 

(3.51) \d^f{t,v)\ < 

for large \t\ and small ry. We proceed by splitting the integration in (|3.48l) and (I3.49P into the 
regions T>buik = [0,1] x {r : \E + r/dwl < 2} and Dout = [ 0,1] x {r : 2 < |H + r/dwl < 4. 
Starting with region Hbuik; the variance bound of Proposition IB. 41 is applicable and we see that 
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K\VN{t + uniformly in t. Combined with (13.511) we see that the integrands of 

(I3.48P and (|3.49l) restricted to Hbuik are dominated by an integrable function and the proof proceeds 
as in the compactly supported case of Theorem 13.21 

It remains to bound the contribution to the integrals on the domain Hout- Here we can exploit 
the decay of the test function / and show that the inner limsup over N will already be zero in 
(j3.48l) and ()3.49p . Therefore it suffices to take H = 0 and K = 1. Then the variance bound of 
Proposition 13.11 applies, yielding 


E(|Hjv(z)|i+^)|5%(t,r?)|'+'^ dtdr, 


1+5 


(3.52) 


J Z)out 


V 


(-l-e)(l+5)+cra|^|(-l-/3)(l-o-)(l+5) 


dtdr] 


< cd'^ 


' Dout 

re(l+5) —(1+;3)(1—cr)(l+(5)+l 


where we choose e and 6 so small that the integral over r] is finite while the integral over t goes 
to zero as N ^ oo. Indeed, if e < aa/3 and 6 < aa/3 then (—1 — e)(l + d) + aa > —1 and the 
r] integral is finite. In the integral over t, we choose a < (5/(1 + J), e < f3/{l + 6) and deduce 
that e(l + (5) — (1 + /?)(! — (t)( 1 + (5) + 1 < 0. To make the bounds work simultaneously we take 
e < min{(TQ:/3,/3/(l + 5)}. We conclude that the limit of (I3.52p is zero. To prove Proposition 13.61 
notice that the integrand is supported on Hout with the same regularity conditions on /. Hence an 
identical calculation to that given in (13.521) shows that (j3.46p converges to zero in as N ^ oo. 
This completes the proof of Propositions 13.71 and 13.61 Consequently, by means of IA.21 this also 
completes the proof of Theorem 13.41 □ 


Corollary 3.8. The sequence of stochastic processes Vm{z) with z G 11 is tight in the space of 
continuous functions on any N-independent rectangle in the upper half-plane H. 


Proof. It suffices to verify the Arzela-Ascoli criterion: 

(3.53) nVN{zi)-VN{z 2 )\^ <C\zi-Z 2 


for zi and Z 2 in some iV-independent rectangle in H and C a constant depending only on the 
vertices of the rectangle (i.e. not on N). To prove this, note that 


(3.54) 


1 1 _ {ui - U 2 ) + i{vi - V 2 ) 

X — ui — ivi X — U 2 — iv 2 {x — ui — ivi){x — U 2 — iv 2 ) 


implies 


(3.55) nVNizi) - VNiz2)\^ = - Z2pE|A^«^°(h)|2 


where h{x) = ((x — ui — ivi)(x — U 2 — iv 2 ))~^ is a smooth function with decay h{x) = 0(|x|“^) as 
|x| —>■ 00 . Hence the techniques of the present subsection are applicable with a = /3 = 1. Indeed, 
after replacing h with a smooth cut-off hjsf as in Lemma 13.5p an application of formula (ll.20p 
followed by Cauchy-Schwarz leads to E|X)0®®°(u)p < where 

1 rl /•oo 

(3.56) 4 :=-/ / ^E| VAr(r + ir])\‘^ r])\ dr drj 

^ Jo J—00 

Following the proof of Proposition 13.71 we easily deduce that 1^ is uniformly bounded in N for fixed 

Ul,Vl,U2,V2 □ 
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Appendix A. Weak Convergence results 

Theorem A.l (Martingale Central Limit Theorem). [HI Theorem 35.12] Let M^^n, ^ ^ k < N, 
N > 1 be a sequence of zero-mean, square-integrable martingales adapted to the filtration iFk,N! 
and let IFo^n denote the trivial a-field. Let denote the martingale difference sequence = 
^k,N — where 1 < A: < A. Suppose the following conditions hold 

N 

(A.l) for all e > 0 0 in probability, 

k=l 

N 

(A.2) at] —)• CT^ in probability, 

k=l 

where Then Mm,n converges in distribution to a Gaussian random variable with mean 0 and 
variance 

Theorem A.2. [THl Theorem 1; Lemma 1] Let ^ be a space of measurable functions 4>{z,w) : 
H X C ^ C such that (j){z, •) is continuous for all z G H. Let M be a set of measurable functions 
a: : El —>■ C such that 

(A.3) J \(l){z, x{z))\dz < Qo, </> G <h, 

and define the functional 

(A.4) = f (j){z,iN{z))dz. 

Jh 

Suppose {^ 7 v : A" G Nq} is a sequence of stochastic processes : El — >■ C, with paths in M . If the 

finite dimensional distributions of converge weakly to those o/^o Lebesgue almost everywhere in 
El and for all e > 0 and (/> G <h.' 

(A.5) lim limsuppf / {\(l){z, {z))\ — K)~^ dz > e] =0, 

K^oo N^oo \Jm J 

(A.6) inf lim sup P ( / \4>{z,f,Niz))\dz>€] =0, 

BcH,A(S)<oo N^oo \Jm-B J 

where A is the Lebesgue measure on H, then for <f>i,... ,4>k ^ we have that 
(A.7) - ■ ■ ,UkiCN)) => (4i(6),---,4fc(Co))- 


Appendix B. Concentration inequalities and bounds on the resolvent 

In this section we record some important bounds required in Section [2l We remind the reader 
the notation used in that Section 

(B.l) 5^"(z) := hlGkizThk - A-i Tr(Gfc(zi)") 

where h^ denotes the kfi^ column of TL with k^^ element removed. The matrix TLk is defined as the 
matrix TL with the k^^ row and column erased and Gk{z) := {TLk — z)~^ is the resolvent of TLk- 
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Lemma B.l. Let z = E + 
constants ci, C 2 snch that 

(B.2) 

(B.3) 


Then for fixed rj > 0 and r G M and all g > 1, there are positive 


z + ^Tv{Gk{z)) 


for all k, N. 

Proof. Ineonalitv (IB.2D can be found in Prooosition 3.2 of m as a consequence of the Hanson- 
Wright large deviation ineqnality m- Such concentration ineqnalities appear in several other works 
concerning the local semi-circle law for Wigner matrices, see e.g. [291 [281 [3Q1 El]. To prove dESj), 
set bk := {z + N~^Tv{Gk))~^ = (z + s{z) + where = A'“^Tr(Gfc) — s{z). Now we nse the 

elementary identity 


(B.4) 


bk ■= 


1 


z -|- s(z) -|- Ak 


E 


At 


^ {z + s{z)y+^ 


+ 




{z + siz))Piz + N-^Tr{Gk)) 


Then it is known that \z + s{z)\ ^ < 1 uniformly on the upper-half plane z E (see Eq. 8.1.19 
in 0). Furthermore, we have the trivial bound \z-\-N ^Tr(GA;)| ^ < {jj/dN + N ^ImTr(GA;)) ^ < 
dN/v- This gives us 


(B.5) 


i6fci <X;iAfcr+^iAfcri 


On the other hand, by rigidity arguments we know that E|Afc|'? < C{Ng/d^) so for example 
E|6fc| <1-1- 0{{Nr]/d]sf)~^) + 0{{d]\f/rj){Nr]/d]sf)~'^~^). This last error term can be made o(l) by 
choosing p large enough (setting djsj = one finds the condition p-|-l>Q;/(l — a)) □ 

Lemma B.2. Let z = E + ■ For all g > 1 with fixed r/ > 0 and r G M, there are positive 

constants c, G such that 


(B.6) 

(B.7) 


E\d],^5^^\zW <G[E] 


E 


^TA{Gl{z)) 


dN J 


< cmax 

\ «Af \ «Af 


rj f Nr] 


-i\ 


for all k, N. 

Proof. Since Gk{z) is an analytic function of z in the upper half plane and £Gk{z) = Gl{z), we 

write 

(B.8) 

= dp^ (^hlGkhk - ^ Tr(Gfc)^ ^ ^ £ {z (^hlGk{w)hk - ^ Tv{Gk{w))^ dw 
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where Sz is a small circle of radius ri/2dj\j around the point z. Note that Sz is of distance 0{r]/dN) 
from the real axis and that \z — w\ = r]/2dN- By Holder’s inequality and Lemma iB.ll we have 


(B.9) 

(B.IO) 

(B.ll) 

(B.12) 




< 


< C 


Nr}\ 

- 9/2 

An J 



-q/2 

An ) 



1 


(27rdAr)'? 


Is, \z-w\ 


;\dw\ 


r] 


As for Tr(G|) we similarly have 


(B.13) ^TiiGl) = ^fTrG, = ^s\z) + 
djsi djsi dz djsi 


Is 


2'KidN J Sz 


{N ^TTGk{w) — m{w))dw 


where Sz is a small circle of radius r]/2d]\f with center 2 : = - 2 (s). Then using the known rigidity 
estimates we get E|A^“^ Tr(Gfc(rc)) — s{z)\'^ < {Nr]/d]^)~‘i which we use to estimate (|B.13ll . Then 


E 






q-i 

E 


{N ^TiGkiw) — m{w))dw 


2ti dN Jsz {w - zf 




< cmax 




as required. □ 

Lemma B.3. The identity ()2.27p holds. 


Proof. We omit the explicit z-dependence on the resolvent and make use of the quantities given in 
5^” in (IB.ip and two further quantities 


(B.14) 


bk ■= 


1 


2 + AT-i Tr Gk ’ 


Gkk = 


'^kk ^ h.Gkhk 
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Using the identity z + N ^ Tr = (—(5^^ + Hkk — we obtain 


(B.15) 

(B.16) 

(B.17) 

(B.18) 

(B.19) 

(B.20) 


l + hjGlhk l + N-^TiGl 
Ukk-z-hlGlhk -z-N-^TiGk 

(I + hjGlhkmkk - G^^il + hjGlhk) ^ G^^{1 + N-^TtGI) 

G^fb^^ 

Gfcfc(l +— 5^^) ck,2r 
- -pi- l^k 

= -blil + hlGlkk^kk - 4') + ^i(i + hiGlhk){nkk - 4')' - ^''/bk 
= -blil + N-^ TiGDiTikk - 5^') - 
- blbfinkk - 4 ') + + hlGlhumkk - 


where to obtain (jB.lSp we expanded using the simple identity Gkk = —bk — bkGkkib^^ — Tikk)- 
Then to obtain ()B.19p from (|B.18p we used that h\Glhk = TrG| + 5^^. The two terms in 
(1B.19D combine as an exact derivative 

(B.21) §^bkinkk - 5^') = -blil + TvGDinkk - 4') - b^/bk 

while the remaining two terms in (jB.20jl combine to give the error term ek^N{z) in (I2.27p . Finally, 
we conclude the proof of the identity by applying (E^ — Efc_i) to (IB.15p . noting that the second 
term vanishes. □ 


Proposition B.4. Fix t) > 0. Then under the same conditions as Theorem \1.^ we have that there 
exists a positive constants Nq, Mq, G, cq, and ci = ci(G, c) such that 

(B.22) E|VAr(z)|2 = Var|d]^iTrG(E + z/d7v)| < ciP^, 

for all N > Nq such that Nrj/djy > Mq, rj/dN < and | i7 + t / d7v | < 2 + 7 // dAr . 

Proof. By the triangle inequality, we can write 
(B.23) 

E|Uv(z)|2 = E\iN/dN)isNiE + z/dn) - Esn{E + z/dN))\'^ 

< 2E| iN/d]\f) (sjv(E' + z/dj\f) — siE + z/d]\f))\‘^ + 2|(iV/d^v )E{s]\fiE + z/d]\f) — s(i7 + z/d^v))P 

< 4E|(Ar/djv)(siv(E + z/dN) - siE + z/djv))P 

where the last line follows from Jensen’s inequality. The latter expectation is by definition the 
integral 

/•OO 

(B.24) / 2MP(|sAr(i7 + z/dAf) — s(E + z/dAr)| > udAr/AI) du 

Jo 

= 7y-2 277P + z/dN) - siE + z/d^v)! > dK 

It suffices to bound the contribution to the above integral when 77 > 1. We apply Theorem 11.31 
with e.g. g = 3 to obtain a convergent estimate. □ 
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Appendix C. Function extension 


The following Lemma is useful for proving results about the linear statistics o of a random 
matrix via analogous results for the resolvent. It can be found in Proof of Lemma 5.5.5 of the book 


Lemma C.l. Let / : M —>■ M be a compactly supported function with first derivative continuous 
(denote it by / E C'^(M)) and consider an extension Tj : —)• C that inherits the same regularity 
and such that = f{x) for all x and Im('I'j(x, 0)) = 0. Further assume that d'i>f{x,y) = 

0{y) as y ^ 0. Then one has 


(C.l) 
where d : 


f{X) + iH[f] = 


vr 


n 


d'i>f{x,y) 
-oo>'-x-iy 


dx dy, 


= ^ + and H : L^(M) —L^(M) is the Hilbert Transform of /: 


(C.2) 


H[f] 



dt 


fit) _ 

X — t 


Remark C.2. Our particular choice of Tj in this paper will be the following (see my- 
(C.3) '^fix,y) = (fix) +i{f{x + y) - f{x))).J{y), 

where J{y) is a smooth function of compact support, equal to 1 in a neighborhood of 0 and equal 
to 0 if y > 1. 


Proof. We make use of the the substitution x ^ x + \ and compute the real part (denoting 
'it f = u + iv) 


(C.4) 

(C.5) 

(C.6) 


1 r°° r 

-Rey dy J 


°° d^f{x + A,y) 

-oo x + iy 


dx 


=-lim 

TT 


dx dy 


X f du{x, y) 
dx 

du{x,y) 


x2 _|_ y2 


X 


x"^ + 


dx 


dv{x,y) \ y / du{x,y) dv{x,y) \ 

dy ) x"^ + y‘^ \ dy dx / ’ 

dvix,y) \ y f du{x,y) dv{x,y) 

dy j + y2 V dy dx 



where = {e < + y^ < R} and R is sufficiently large (by the compact support hypothesis). 

Changing to polar coordinates in the latter integral, a simple computation shows that the integral 
transforms as 


(C.7) 

(C.8) 

(C.9) 


-lim 

vr 


pn pR 

im / d6 I dr 

Jo Jt 


9u(rcos(0), rsin(0)) d i;(rcos(6*), rsin(0)) 


dr 


de 


- — lim ® 


/ u(rcos(0),rsin(0)) , 

u{r cos(0), r sm(0)) dO H- - dr 


lim[ [ dO u{R cos(6), R sm(9)) — [ d6 ule cos(9), e sm(9)) + [ 

TT ^^0 \ Jo Jo Je- 


v{r, 0) 

<|r|<R X 


dr 


where in the second line we applied Green’s theorem to reduce the double integral to an integral 
over the boundary of the positively oriented rectangle Re,R with vertices (e, 0), {R, 0), {R, vr), (e, vr). 
Now choosing R large enough and using the assumption Im('I'j(r + A, 0)) = u(r, 0) = 0 we see 
that the above limit is equal to rt(0,0) = 'I'j(A,0) = /(A). The proof for the imaginary part is 
similar. □ 
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